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ABSTRACT 


The  underwater  spark  discharge  is  a  means  of  converting  electrical 
to  acoustic  energy.  In  this  process,  the  medium  is  a  fundamental  part  of 
the  transduction  process.  An  underwater  spark  discharge  is  formed  by 
switching  an  energy  storage  capacitor  at  a  high  voltage  to  a  pair  of  sub¬ 
merged  electrodes.  An  ionized  plasma  at  high  temperature  and  pressure 
is  formed  between  the  electrodes.  The  pressure  within  the  plasma  is 
transmitted  to  the  surrounding  water  as  a  wave  of  compression  similar 
to  the  shock  wave  created  by  an  underwater  explosion. 

An  analysis  of  the  electro-acoustic  properties  of  the  underwater 
spark  is  presented.  The  analysis  is  based  upon  expressing  the  macro¬ 
scopic  properties  of  the  gaseous  plasma  in  terms  of  equations  for  the  con¬ 
servation  of  mass  and  momentum.  The  non- compressive  approximation 
is  used  together  with  the  assumption  of  adiabatic  expansion,  to  derive  and 
solve  an  equation  of  motion  for  the  plasma  radius  as  a  function  of  time. 
From  the  radius -time  relationship,  formulas  for  the  plasma  pressures 
and  spark  channel  wall  velocity  are  derived. 

The  sonic  energy  radiated  from  the  underwater  spark  discharge  is 
discussed,  and  the  directional  properties  of  the  radiated  pressure  are 
derived. 

Properties  of  the  underwater  spark  discharge  plasma  as  a  circuit 
element  are  investigated.  An  electrical  equivalent  for  the  discharge 
circuit  is  developed  from  a  knowledge  of  the  variation  of  channel  area 
with  time.  The  time  varying  plasma  resistance  is  used  in  the  differential 


equation  of  the  circuit,  and  solutions  are  obtained  for  the  current  and 
voltage  waveforms. 

The  theoretical  results  are  compared  with  experimental  for  a  fresh 
water  underwater  spark  discharge  having  the  following  parameters, 
electrode  voltage,  v  =  30  kilovolts,  storage  capacity,  C  =  .  1  microfarads, 
and  gap  separation,  L  =  4  cm.  Considering  the  nature  of  the  approxi¬ 
mations  made  in  the  development  of  the  theory,  good  agreement  with 
experimental  results  are  obtained.  In  particular,  the  theoretical  de¬ 
pendence  of  the  spark  channel  radius  upon  time  is  found  to  agree  with  the 
experimental  results  of  several  investigators. 

The  theory  developed  has  resulted  in  an  improved  understanding  of 
the  processes  encountered  in  an  underwater  spark  discharge.  Quantities 
which  are  extremely  difficult  to  measure  in  the  laboratory  may  now  be 
predicted.  For  example,  it  is  found  that  for  the  spark  discharge  para¬ 
meters  stated  above,  the  spark  channel  formed  has  an  initial  radius  of 
about  .4  millimeters  and  an  internal  pressure  of  IZ,  000  atmospheres. 
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I.  INTRODUCTION 

An  underwater  spark  discharge  is  produced  when  the  energy 
from  a  high  potential  source  is  rapidly  transferred  to  a  set  of  submerged 
electrodes.  Upon  the  breakdown  of  the  water  gap,  a  plasma  of  ionized  gas 
forms  between  the  electrodes  resulting  in  a  rapid  discharge  of  the  energy 
source.  Since  the  inertia  of  the  surrounding  water  exerts  an  almost  rigid 
opposition  to  the  expansion  of  the  spark  channel  formed,  extremely  high 
pressures  are  developed.  Investigators  have  reported  channel  pressures 
in  excess  of  10,  000  atmospheres  and  temperatures  greater  than  40,  000 °K. 

A  schematic  diagram  of  the  type  of  electrical  circuit  and 
electrode  assembly  used  to  produce  an  underwater  spark  discharge  is  shown 
in  figure  1.  When  the  switch  tube  is  fired,  the  high  voltage  is  transferred 
to  the  electrodes.  The  high  intensity  electric  field  developed  across  the 
gap  causes  the  ionization  of  the  water  molecules  and  the  formation  of  the 
gaseous  plasma. 

Edward  A.  Martin^  reports  upon  an  extensive  investigation  into 
the  properties  of  the  discharge  plasma.  A  Kerr  Cell  photograph  made  by 
Martin  of  a  25  kilovolt  spark  discharge  plasma  at  4.  1  microseconds  is 
shown  in  figure  2. 

Through  electrical  measurements  made  upon  the  discharge 
circuit  and  Kerr  Cell  photographs  of  the  plasma,  Martin  was  able  to  cal¬ 
culate  the  energy  distribution  among  the  circuit  components  and  plasma  at 
different  times  during  the  discharge.  The  plasma  energy  calculations  were 
then  compared  to  the  energy  which  the  constituent  particles  of  the  spark 
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channel  would  store  at  the  experimentally  determined  temperature  and 
pressure. 
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The  electro-acoustic  efficiency  of  the  underwater  spark  dis¬ 
charge  has  been  investigated  by  two  Soviet  scientists,  N.  A.  Frolov  and 
2 

D.  P.  Roi  .  Their  results  have  indicated  that  electro-acoustic  ef¬ 
ficiencies  of  up  to  32  percent  can  be  achieved  with  spark  gap  lengths  of 
5  cm  at  voltages  of  30  kilovolts. 

In  many  respects,  the  underwater  spark  discharge  possesses 

the  same  hydrodynamic  and  acoustic  properties  as  a  small  scale  under- 
3  4 

water  explosion  ’  .  In  an  underwater  explosion,  a  gas  globe  consisting 

of  the  detonation  products  is  formed  at  high  temperatures  and  pressures. 

The  inertia  of  the  water  and  the  elastic  properties  of  the  gas  and  water 

together  provide  the  necessary  conditions  for  an  oscillating  system  and 

the  gas  globe  undergoes  repeated  cycles  of  expansion  and  contraction. 

The  analogy  between  the  gas  globe  formed  by  an  underwater  explosion,  and 

the  underwater  spark  discharge  plasma  has  been  discussed  in  a  published 

5 

work  by  the  author  ,  and  will  be  further  demonstrated  in  this  dissertation. 
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II.  FORMATION  OF  THE  GASEOUS  PLASMA 

Following  the  application  of  a  high  potential  to  the  submerged 
electrodes,  the  sequence  of  events  as  reported  by  Martin^  is  as  follows; 

A  formative  period  of  the  order  of  a  microsecond  exists,  during 
which  branched  streamers^  of  ionized  molecules  advance  from  each 
electrode  until  one  streamer  makes  contact  with  the  other  electrode.  The 
discharge  then  takes  place  through  this  streamer  causing  the  others  to 
conduct  negligible  current. 

The  duration  of  the  formative  period  may  be  observed  in  the 
current  and  voltage  oscillograms  taken  from  Frolov's  and  Roi's  paper  and 
reproduced  in  figure  3.  For  example,  with  a  gap  length  of  5  cm,  ap¬ 
proximately  2-1/4  microseconds  following  the  application  of  30  kilovolts 
to  the  gap,  breakdown  occurs.  The  potential  across  the  gap  drops  abrupt¬ 
ly  and  the  rate  of  rise  of  current  increases.  As  the  gap  length  is  reduced, 
both  the  formative  time  and  the  time  between  breakdown  and  maximum  cur¬ 
rent  are  reduced.  With  a  reduction  in  the  gap  length,  the  initial  drop  in 
the  voltage  becomes  even  greater. 

At  the  instant  of  breakdown,  a  gaseous  plasma  at  a  high  pres¬ 
sure  and  temperature  is  formed  between  the  electrodes.  It  is  evident  that 
these  conditions  cannot  persist  without  affecting  the  surrounding  medium. 
Due  to  the  continuity  of  pressure  through  the  plasma-water  boundary,  a 
high  pressure  must  also  exist  in  the  water  a.djacent  to  the  spark  channel. 
This  abrupt  increase  in  pressure  when  the  plasma  is  formed  is  propagated 
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outward  as  a  wave  of  compression  in  the  water. 

The  pressure  magnitudes  are  such  that  acoustic  theory  does 

not  apply  to  the  propagation  of  such  a  wave.  This  type  of  steep-fronted 

wave  whose  pressure  amplitude  is  sufficiently  great  to  change  the  density 

7 

of  the  medium  is  commonly  called  a  "shock"  wave  .  A  detailed  treatment 
of  the  propagation  of  the  cylindrical  shock  wave  generated  by  the  under¬ 
water  spark  discharge  is  extremely  complex  and  will  not  be  attempted  in 

this  dissertation.  The  interested  reader  may  be  referred  to  several  of 

^  8,  9 

the  references 


III.  ACOUSTIC  PROPERTIES  OF  THE 
UNDERWATER  SPARK  DISCHARGE 


1.  Basic  Relationships 

The  acoustic  properties  of  the  underwater  spark  discharge  may  be  de¬ 
rived  from  the  equation  of  motion  of  the  gaseous  plasma.  Solution  of  the 
equation  of  motion  will  yield  the  motion  of  the  plasma  boundaries  with 
time.  A  knowledge  of  this  motion,  together  with  the  pressure- volmne 
relationships,  enable  the  pressure  within  and  on  the  surface  of  the  plasma 
to  be  calculated. 

An  adequate  description  of  the  motion  of  the  plasma  may  be  developed 
if  it  is  assumed  that  the  non-linear  nature  of  the  propagation  does  not  af¬ 
fect  the  energy  relationships  which  govern  the  differential  equations  of 
motion.  With  this  approximation,  the  energy  relationships  may  be  de¬ 
rived  on  the  basis  of  linear  acoustic  theory.  This  technique  has  been 
used  successfully  by  several  authors  in  the  formulation  of  the  equations 
of  motion  for  the  gas  bubble  formed  by  an  underwater  explosion,  where 
pressures  exist  which  are  an  order  of  magnitude  greater  than  the  under- 

^  1  J-  V.  10,  11 

water  spark  discharge 

We  will  use  the  equations  for  the  conservation  of  mass  and  momentum 

for  an  ideal  fluid  medium.  These  equations  may  be.  obtained  from  any 
12  13 

standard  text  ’  on  acoustics  or  hydrodynamics  and  are  summarized 


below  for  convenience. 
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Conservation  of  Mass 


Let  be  the  density  of  the  fluid,  and  v,  the  vector  velocity  of  a 
point  moving  with  the  fluid;  then  the  equation  for  the  conservation  of  mass 
is; 


div  (/O  v) 


d  t 


Conservation  of  Momentum 

If  p  is  a  scalar  pressure  distribution  which  may  be  a  function  of 
time  and  the  spatial  coordinates;  then  Newton's  second  law  requires  that, 


dv 

(2)  grad  p  =  - /O 


A  considerable  simplification  of  these  equations  results  when  the 
non- compressive  approximation  is  made.  Physically,  this  approximation 
is  the  assumption  that  the  pressure  variations  are  sufficiently  small  as  to 
cause  a  negligible  change  in  the  density  of  the  medium.  The  non-linear 
effects  due  to  finite  compressibility  will  be  discussed  further  when  the 
propagation  of  the  shock  wave  due  to  the  underwater  spark  discharge  is 
discussed  in  Section  3.  12. 


Using  the  non-compressive  approximation,  remains  unchanged 
with  time  ,  and  we  may  modify  (1)  accordingly, 


(3) 


div  V  =  0 


In  acoustics,  just  as  in  electromagnetics,  (3)  implies  the  existence 
of  a  scalar  potential  ^  whose  negative  gradient  is  equal  to  v. 

(4) 

The  scalar  velocity  potential  must  be  a  solution  to  Laplace's 
equation  which  may  be  obtained  by  substituting  (4)  in  (3); 

(5)  =  0 

Examining  (2)  we  see  that  if  we  substitute  (4)  for  v  we  obtain; 

VP=-^  -V^) 

or, 

(6)  P  =  /? 

Which  is  the  relationship  between  pressure  and  velocity  potential  in 


an  incompressible  fluid. 


2.  Motion  of  the  Gaseous  Plasma 


Before  the  equation  of  motion  of  the  gaseous  plasma  may  be  form¬ 
ulated,  it  is  necessary  to  make  several  other  assumptions; 

It  is  assumed  that  at  the  instant  of  breakdown  of  the  gap,  a  gaseous 
plasma  having  an  instantaneous  temperature  and  internal  pressure  p^ 
is  formed.  Furthermore,  the  shape  of  the  plasma  is  such  that  it  may  be 
considered  a  right  circular  cylinder  of  length  equal  to  the  gap  length  L, 
and  initial  radius,  a^.  A  uniform  energy  density  throughout  the  plasma 
is  also  assumed. 

A  diagram  of  the  plasma  and  the  coordinate  system  which  will  be  em 
ployed  is  shown  in  figure  4. 

In  addition,  the  assumption  is  made  that  from  the  instant  of  break¬ 
down,  some  fraction  of  the  original  capacitor  energy  is  available  for 
plasma  motion.  Actually,  this  energy  is  supplied  at  the  start  of  current 
flow.  However,  if  the  plasma  motion  is  slight  during  the  interval,  there 
is  little  error  in  treating  the  energy  as  available  from  the  beginning  of 
the  interval.  As  a  first  approximation,  the  radiation  of  acoustic  energy 
will  be  neglected,  as  well  as  dissipation  of  energy  by  thermal  conduction 
at  the  plasma- water  boundary. 


The  kinetic  energy  of  motion  is  the  energy  associated  with  the  water 
by  reason  of  its  motion.  If  spherical  instead  of  cylindrical  symmetry 
existed,  this  energy  could  be  derived  quite  simply.  We  shall  then  modify 
the  equations  to  take  cylindrical  symmetry  into  account. 


Let  u  =  dr/dt  be  the  radial  velocity  of  the  fluid  at  a  point  r.  Then,  by 
virtue  of  equation  (3) 


Upon  performing  the  partial  differentiation, 

2ru  +  r  ^ 
or, 

(8)  du  _  2u 
dr  ~  r 

Applying  the  boundary  condition  that  at  the  plasma  wall 
r  =  a,  u^  =  da/dt  and  solving  for  u, 

(9)  u  =  — -  d  a/dt  r  >  a 

Hence,  in  the  water,  the  particle  velocity  falls  off  as  the  square  of 

the  radial  distance.  Therefore,  for  the  spherical  case,  the  kinetic  energy 
of  flow  is , 


14 


K.  E. 


^  (4'7fr^dr) 


da 


K.  E.  g  =  2  Tf/o  a' 


(-*-) 


Now,  consider  a  cylindrical  plasma  of  radius  a  and  length  L.  In  the 
cylindrical  region  a  <  r  <  L/2,  the  plasma  approximates  an  infinite 
cylinder,  therefore  the  velocity  potential  jp  falls  off  approximately  as 

— — “j—  *  So  that  by  equation  (4),  the  velocity  will  fall  off  approximately 


Employing  the  boundary  condition,  we  find  for  this  region  the  fol¬ 
lowing  approximate  expression  for  the  particle  velocity; 


a  r  <r  L/2 


Therefore,  the  kinetic  energy  for  the  region  under  consideration  is, 


K.  E. 


nr/oa  L 


i^y 
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Since  the  plasma  radius  is  assumed  to  be  much  less  than  L/2,  the* 
value  of  the  K.  E.  at  L/2  will  be  neglected,  as  well  as  any  contributions 
for  r  >  L/2.  Therefore,  the  kinetic  energy  of  the  water  is; 


(11) 


K.  E. 
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4.  Potential  Energy 


The  potential  energy  of  the  spark  channel  is  the  work  done  by  the 
channel  in  expanding  from  a  radius  a  to  infinity. 

.e/> 

(  p  -  P)dv  = 


(12) 


>.  E.  =  ( { 

J 


(p  -  P  )  Zir  aLda 


a  a 

P  is  the  hydrostatic  pressure  at  the  depth  of  the  discharge,  and  p  is 
the  external  plasma  pressure.  A  further  physical  interpretation  of  the 
potential  energy  as  well  as  the  external  plasma  pressure  will  be  given 


later  in  Section  3.  7. 
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5.  Equation  of  Motion 


We  have  assunaed  that,  from  the  instant  of  breakdown,  the  plasma 
possesses  an  energy  Y  which  remains  constant.  On  the  basis  of  this 
assumption,  an  equation  of  motion  may  be  written  which  expresses  the 
conservation  of  kinetic  and  potential  energy.  From  (11)  and  (12); 

(13) 

a 


oo 


(p  -  P)  ZTfaLda  =  Y 


Using  the  fact  that  at  the  beginning  of  the  motion  when  a  =  a  , 


da/ dt  =  0; 


(14) 


.  CO 


Y  =  /  (p  -  P)  Zra  JLda 

a 


Y  has  the  physical  significance  of  being  the  fraction  of  the  initial  capacitor 
energy  which  is  available  for  plasma  motion. 

Substituting  Equation  (14)  in  Equation  (13)  we  obtain  the  following 
equation  of  motion; 


(15) 


2 

/O 


(p  -  P)  a  da 


a 


This  equation  expresses  the  relationship  between  the  plasma  radius, 
wall  velocity,  and  wall  pressure. 

Equation  (15)  may  be  integrated  provided  the  dependence  of  p  on  a  is 
known.  If  it  is  assumed  that  in  the  expansion  of  the  plasma,  no  heat  is 

transferred  to  the  surrounding  water,  an  adiabatic  law  is  followed. 


6.  Adiabatic  Expansion 

We  assume  that  the  underwater  spark  discharge  plasma  behaves  as  an 

ideal  gas.  This  is  not  strictly  true  due  to  the  production  of  long  range 

interparticle  forces  by  the  high  degree  of  ionization.  However,  we  will 

use  the  ideal  gas  law  as  a  first  approximation,  as  Martin  does  in  his  cal- 

14 

culations.  When  an  ideal  gas  undergoes  a  quasistatic  adiabatic  process; 


Cpd^  -  vdp  =  0 


where,  c  =  specific  heat  at  constant  volume 
Cp  =  specific  heat  at  constant  pressure 
d^  =  differential  temperature  change 


Dividing  Equations  (16) 

dv 
v 

c 

or,  regarding  )["  =  — —  as  constant  and  integrating; 

v 

y  y 

(17)  pv  =  p  V  =  constant 

^  ^o  o 

The  ratio  of  the  specific  heats  ^  maybe  obtained,  provided  we  as¬ 
sume  that  the  plasma  formed  by  the  underwater  spark  discharge  behaves 
as  a  fully  ionized  gas.  Plasmas  of  5  percent  ionization  are  generally  con¬ 
sidered  completely  ionized,  meaning  that  the  long  range  interparticle 
Coulomb  forces  play  a  dominant  rolo  in  determining  the  plasma  properties. 


Martin  reports  that  the  underwater  spark  discharge  is  approximately  30 

15  .... 

percent  ionized  ;  so  that  the  assumption  of  "complete"  ionization  is  a 

reasonable  one. 


For  a  fully  ionized  gas,  Spitzer^^  gives  the  following  formula  for  ^  ; 


(18) 


2  +  m 
m 


v,'here  the  quantity  m  is  the  number  of  external  degrees  of  freedom. 

Since  a  heavy  current  flows,  particle  motion  exists  in  a  direction  parallel 
to  the  electrodes.  In  addition,  the  cylindrical  plasma  is  expanding  radially, 
so  that  the  particles  have  a  radial  motion.  This  argument  would  seem  to 
imply  the  existence  of  two  external  degrees  of  freedom  resulting  in 
y  =  2.  If  the  particle  velocities  were  randomized,  due  to  collisions, 
y  would  equal  5/3  which  is  sufficiently  close  to  ^  =  2  to  warrant  the  use 
of  the  value  y  =  2  in  subsequent  calculations. 


The  solution  obtained  by  the  adiabatic  assumption  expresses  the  kinetic* 
pressure  of  the  plasma  as  a  function  of  the  plasma  volume.  In  the  solution 
of  the  equation  of  motion,  the  external  pressure  (the  pressure  which  ex¬ 
ists  in  the  water  adjacent  to  the  spark  channel)  must  be  used.  The  as¬ 
sumption  will  be  made  that  the  kinetic  and  external  pressures  are  equal. 

This  is  not  strictly  true  due  to  deviations  from  the  ideal  gas  law.  More 
will  be  said  about  the  plasma  pressures  in  Section  3.  8. 


*The  pressure  used  in  the  ideal  gas  law,  p,  .  v  =  MB.T 

®  ^kinetic 


7.  Solution  of  the  Equation  of  Motion 


Using  the  adiabatic  approximation  for  the  external  pressure,  the  plasma 
radius  as  a  function  of  time  will  now  be  obtained.  Substitution  of  the  de¬ 
pendence  of  the  external  pressure  p  upon  the  plasma  volume  from  (17)  into 
(15)  results  in  the  following  integrodifferential  equation  for  the  radius  of 
the  underwater  spark  discharge  plasma  as  a  function  of  time; 


(-4^)  f  bo  (^)  -s] 


We  may  express  the  variable  a  in  terms  of  the  non-dimensional  vari¬ 
able  x  =  The  equation  becomes; 


p^  X  ^  "  P  X  dx 


Performing  the  integration,  we  obtain; 


C  1  -  x'^  +  o<(  1  -  x^)  J 


Where  oC  =  P/p^>  is  generally  <<  1  during  the  early  stages  of 


the  motion. 


Separating  the  variables  in  (21)  and  solving  for  t, 


t  =  a 


x^  dx 


o  /  p 

V 


Px^  -1  +  (1  -  x‘ 


Since  the  early  stages  of  the  motion  when  1  are  of  interest,  the 

term  involving  ©<'  may  be  neglected  and  the  integration  of  (22)*  results  in; 


t  =  X  .  /  - 1  +  cosh 


'  -J 


whe  r  e , 


J  “Po 


This  is  the  equation  for  the  variation  of  the  plasma  radius  with  time. 
An  approximation  to  this  expression  for  a/a.^  >  1.  2  is; 


t  «  ^x^ 


a  ^  ao  (t/'t  ) 


Equation  (23)  and  the  approximation  (24)  are  plotted  in  figures  5  and  6. 


An  alternate  form  of  (25)  may  be  derived.  Using  the  adiabatic 
equation,  we  may  obtain  a  solution  to  (14). 


The  total  energy  Y  is  given  by; 


(p  -  F)  2'7r  aLda 


2  If  aLda 


p  TTa**  L 
^o  o 


^Pierce's  Table  of  Integrals  P.  22  #150 


Substituting  (26)  into  (25),  an  equivalent  expression  for  the  underwater 
spark  discharge  plasma  radius  as  a  function  of  time  is  obtained; 


TT/O  L 


a  »  a 


Equation  (27)  may  be  derived  a  bit  more  directly  by  a  simple  physical 
interpretation  of  the  equation  of  motion.  Let  us  examine  the  equation  for 
the  potential  energy  of  the  plasma  obtained  from  (12); 


P.E.  = 


pdv  +  P  TT  a^L 


The  second  term  on  the  right  represents  the  work  done  by  the  plasma 
against  hydrostatic  pressure.  The  first  term  represents  the  internal 
energy  of  the  plasma  as  a  function  of  its  radius; 


E(a)  = 


pdv  =  p  V 
^  o 


p  v^ 
^o  o 


This  shows  that  the  internal  energy  decreases  as  rapidly  as  a*;  so 

that  for  a  >>  a^,  the  internal  energy  becomes  negligible  compared  to  the 

kinetic  energy.  In  the  same  manner,  if  a  «  a  ,  where  a  is  the  maxi- 

m  m  , 

mum.  radius  of  the  plasma,  since  p  >>  P,  the  hydrostatic  term  may  be 
neglected  in  comparison  to  the  kinetic  energy.  We  are  then  left  with  an 
approximate  equation  of  motion, 


r  j 

a  «  a  «  a 
o  m 
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Integrating  (31); 

Equation  (32)  is  identical  to  (27).  Both  equations  show  that  the  spark 
discharge  plasma  radius  varies  as  the  square  root  of  the  time. 


1/4 


1/2 


a 


<r< 


m 


Spark  Channel  Wall  Velocity 

The  wall  velocity  of  the  expanding  spark  channel  may  be  obtained 
from  (21).  Neglecting  <><■  ,  the  channel  wall  velocity  in  the  early  stages 
of  the  motion  is, 


1/2 


a 

o 


-1  J 


1/2 


The  maximvim  wall 
It  occurs  when  a  =  ^  2  a^ 


velocity  may  be  obtained  by  differentiating  (33). 
and; 


(34) 


da 

"3r 


MAX 


8.  Plasma  Pressure 

In  the  solution  of  the  equation  of  motion,  we  have  assumed  that  the 
underwater  spark  discharge  plasma  obeys  the  ideal  gas  law.  This  law 
concerns  a  given  amount  of  a  gas  at  a  temperature  T.  The  gas  is  assumed 
to  have  no  interparticle  forces  as  well  as  particles  of  zero  volume.  The 
average  pressure  which  is  required  to  constrain  these  particles  to  a  given 
volume  is  designated  as  the  kinetic  pressure  of  the  gas 
ideal  gas  law  applies  to  this  kinetic  pressure  and  states  that; 

(35)  p  V  =  RT 

kinetic 

where  the  gas  constant  R  is  related  to  the  universal  gas  constant  in  a 
mixture  of  j  different  masses  ^  j  by, 

(36)  R  =  Rm^ 

j 

Mj  being  the  molecular  weight  of  the  jth  component. 

The  total  kinetic  pressure  is  the  sum  of  the  partial  pressures  and; 

Pkina.ic  =  NKT 

where, 

N  =  number  of  particles  per  unit  volume 

-Z3  “ 

K  =  Boltzman's  constant  =  1.  38  (10)  joules/  K 

If  the  underwater  spark  discharge  plasma  behaved  as  an  ideal  gas, 
the  external  pressure  necessary  to  constrain  the  gas,  p  .  would 

equal  j However,  the  existence  of  a  high  degree  of  ionization 

produces  long  range  interparticle  forces  resulting  in  additional  pressures. 


If  the  particles  exert  repulsive  forces  on  each  other,  then  the  tend¬ 
ency  to  expand  is  increased.  The  increase  in  the  constraining  pressure  is 

called  p  ,  .  .  If  there  is  also  an  interparticle  attraction,  the  pres- 

^repulsion  ^  ^ 

sure  on  the  constraining  envelope  is  reduced,  the  reduction  being  called 

p  . .  ...  .  It  is  usual  to  call  the  difference  of  attractive  and  repulsive 

^attraction  ^ 

pressures  p.  ,  ,,  which  is  usually  positive. 

^  ^internal  ^ 

In  addition,  if  there  is  a  motion  of  charged  particles  through  the  gas, 
the  resultant  magnetic  field  produces  average  forces  tending  to  force  the 
particles  together.  This  reduces  the  required  envelope  pressure  by  an 
amount  designated  as  More  will  be  said  about  this  pinch  pressure 

when  the  electrical  properties  of  the  discharge  are  discussed  in  Section 
4.  2. 


The  overall  effects  of  these  pressure  contributions  may  be  summar¬ 
ized  in  the  following  equation, 

^kinetic  ^external  ^  ^pinch  ^  ^internal 
17 

Martin  has  found  that  the  internal  pressure  has  a  maximum 
magnitude  of  about  1/4  of  the  external  pressure  which  he  obtained  by 
means  of  the  measured  wall  velocity  and  shock  wave  equations. 

It  is  to  be  noted  that  in  the  evaluation  of  the  equation  of  motion,  the 
external  pressure  was  assumed  to  be  equal  to  the  kinetic  pressure. 
However,  the  kinetic  pressure  is  greater  than  the  external  pressure  due 
to  the  pinch  and  internal  pressure  contributions.  This  increase  will  tend 
to  offset  the  error  due  to  the  noncompressive  approximation  in  which 
lower  external  pressures  than  actually  exist  are  assumed. 
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The  kinetic  pressure  (or  external  pressure)  may  be  evaluated  by 
means  of  the  adiabatic  expansion  equation  (17); 


(39) 


external 


Pk 


inetic 


Using  (23)  for  x(t),  a  plot  of  p/p^  vs  t/^  is  shown  in  figure  7.  The 
similarity  of  this  pressure  curve  to  an  exponential  decay  is  evident.  In 
fact,  an  exponential  approximation  will  be  used  in  the  next  section.  This 
approximation  is  shown  in  figure  7  as  a  dotted  curve  and  may  be  written  as, 


(40) 


Pkinetic^Po 


•EXACT  EXPRESSION 
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9.  The  Directional  Distribution  of  the  Acoustic  Radiation 

The  directional  distribution  of  the  acoustic  radiation  from  an  under¬ 
water  spark  discharge  will  now  be  derived  using  an  acoustic  approximation. 
The  plasma  column  can  be  considered  as  a  line  radiator  having  a  length 
equal  to  the  gap  separation  L.  For  a  point  R  sufficiently  far  away,  lines 
from  different  points  may  be  considered  parallel.  Figure  8  depicts  the 
geometry. 


(42) 


(t,  )  = 


—  oo 


31 


where, 


u  -  y/L 

^  ^  =  L/c  sin  ^ 

and  +1/2 

a  (u)  du  =  1 

-1/2 


In  the  special  case,  where  f  (t)  =  ^(t).  we  denote  the  response  by 


hj^  (t.  ); 


(43) 


^(t  +  a  (u)  du 


—  OA 


Using  the  sifting  property  of  the  delta  function, 


(44)  hj^{t,  4r)  = 


-a  (t/AY)  It  I  <  l^'t'/Zf 
It  I  >/Ar/2/ 


A't' 

o 


Therefore,  the  conclusion  is  reached  that  the  plane  wave  impulse  re¬ 
sponse  hj^  (t,  ^f)  is  the  same  function  of  (t/A't')  as  the  strength  distribution 
is  of  u. 


Using  the  properties  of  convolution,  the  waveform  at  R  due  to  a  line 
radiator  having  an  amplitude  distribution  a  (u)  is, 


(45) 


Pp^  (t, /ITT) 


r 


-a  (-  z/-d  ^  ) 


f  (t  -  z)  dz 


i-oo 


V*  e  will  now  use  equation  (44)  for  an  evaluation  of  the  directional  dis¬ 
tribution  of  the  pressure  due  to  an  underwater  spark  discharge. 


If  we  assume  that  the  strength  of  the  radiated  wave  is  uniform  over  the 
length  of  the  discharge  plasma; 

(46)  a  (u)  =  U  (u  +  1/2)  -  U  (u  -  1/2) 


Inserting  the  above  and  the  exponential  approximation  e 
into  (45),  we  obtain; 

o<:(t  -  z) 


-  oC  t  for  f(t) 


(47) 


Pr  (t,  = 


e 


U  (t  -  z  )  dz 


where, 


=  the  reciprocal  time  constant 


Equation  (47)  is  solved  in  Appendix  A.  The  result  is. 


(48) 


(t,  ) 


s inh  / 2  -oct 

tt.A'i  12  ® 


t 


l/oc.'t'  C  1  -  e-^^e  - 


I2<t  /2 


t  ^  -^'t/2 


This  solution  is  plotted  in  figure  9  as  a  function  of  o<  t  for  «'!_;/ c  =  5 
and  ^  =  0°,  15°,  30°,  60°,  and  90°.  At  9  -  0° ,  the  waveform  is 
identical  to  f  (t)  =  e  since  all  contributions  arrive  simultaneously. 

At  5  =  90°,  since  oc  L/ c  is  large  compared  to  the  time  constant  of  the 
waveform,  the  pressure  at  R  approaches  the  same  form  as  a  (u).  f  (t) 
may  then  be  regarded  as  an  impulse. 


The  peak  value  of  p^  (t,  A't )  occurs  at  t  =  +  ^'^/2  and  is, 

K. 


(49) 


PRpeak(^'"/2.Ar)=  C-1 


-  otA^ 


oiA'C 


Equation  (49)  is  plotted  in  figure  10  as  a  function  of  Q  with  o^L/ c  as  a 
parameter.  Since,  physically,  ®^L/c  represents  the  ratio  of  the  travel 
time  across  the  source  to  the  time  constant  of  the  waveform,  in  limit  as 
o(:L/c-~>0,  the  directional  dependence  vanishes. 


1 0 .  Acoustic  Energy  Flux 


The  acoustic  energy  radiated  by  the  underwater  spark  discharge  will 
be  discussed.  It  was  assumed  earlier  that  there  was  no  radiation  of 
acoustic  energy.  If  the  solution  obtained  is  taken  as  a  first  approximation, 
it  is  possible  to  calculate  what  energy  is  actually  radiated.  First,  the 

•V. 

energy  flux  normal  to  the  spark  channel  will  be  calculated.  If,  in  ad¬ 
dition  to  this  calculation,  the  directional  distribution  of  radiated  energy 
is  obtained,  the  total  energy  flux  radiated  by  the  underwater  spark  dis¬ 
charge  plasma  will  be  known. 

The  instantaneous  energy  flux  per  unit  area  per  unit  time  in  any 
direction  is, 


(50) 


d  0 

3F 


/O  c 


where  /Oc  equals  the  wave  impedance  of  the  medium. 


If.  we  consider  a  direction  normal  to  the  spark  channel,  where 
^  =  0°,  and  assume  first  cylindrical  and  then  spherical  spreading,  the 
pressure  at  a  distance  r  is  approximately. 


(51)  P, 


ext. 


1 


normal 

Substituting  (5  1)  into  (50),  and  separating  variables. 


.orra.l  = 

2  r  ^c 


^a  L  A  ^  =  ^1 


y 

X  =  1 


dt 

■~8~ 

X 


If  dt  is  substituted  from  (21) 
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(53) 


and, 


(54) 


dt  = 


(x^  -1) 


TTZ 


a  j-^ —  dx 


X  =  X 


3/2 


o  L 


normal  ^  2-1/2 


r  /O  c 


dx 

x‘’(x2  - 1) 


X  =  1 


Since  the  major  part  of  the  integral  arises  from  x  near  unity,  little 
error  will  be  introduced  by  letting  Xj  =oo  .  The  integral  may  be  easily 

:}c 

evaluated  as  4/5  .  2/3  =  8/15. 


The  expression  for  p  ,  becomes, 

^  /  normal 


(55) 


8 

normal  “W 


3/2  2  , 

p  a  Li 
_ o 


It  is  interesting  to  check  this  result  with  one  obtained  by  use  of  the 

_t/2  T 

exponential  approximation  for  the  external  pressure,  poCe 


(56) 


normal 


oo 


y 


r 


dt 


'  normal 


)  2  L 
o  o 


2 


e  at; 


p  a  L 

°  °  .  r 

2  /O  c 


Successive  use  of  #205  Pierce’s  Table  of  Integrals,  Ginn  &  Co.  1929 
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The  normal  energy  flux  becomes  , 


Therefore,  an  error  of  about  3%  is  introduced  by  substituting  the  ex 
ponential  approximation  for  the  pressure  in  the  calculation  of  the  energy 
flux  normal  to  the  spark  channel. 


1 1.  Directional  Distribution  of  Radiated  Energy 


The  energy  flux  radiated  in  a  direction  normal  to  the  spark  channel 
has  been  calculated.  The  energy  flux  radiated  in  directions  other  than 
the  normal  may  be  calculated  by  means  of  the  directional  distribution  of 
the  radiated  pressure.  An  energy  flux  directivity  factor  may  be  defined 
as  the  ratio  of  the  integral  of  the  energy  flux  over  a  sphere  to  the  normal 
energy  flux; 


(58) 


'tr 


^(0,^  )sin^dfl  d/ 
^  normal 


If  o<  L/ c  is  sufficiently  small  so  that  ^ 

normal 

=  1  and  the  source  is  essentially  spherical. 

In  terms  of  ^  the  total  radiated  acoustic  energy  W  is, 


then 


(59) 


W 


'H'  Po  ^O  ^ 


D  /  normal 


rrrc 


‘i 


D 


The  acoustic  efficiency  may  be  defined  as  W/Y.  Using  (26)  for 


Y; 


The  result  (6l)  is  significant  in  that,  since  in  the  acoustic  approxi¬ 
mation  there  are  no  energy  losses  with  range,  is  the  fraction  of 

Y  actually  radiated  as  acoustic  energy.  The  energy  measurable  over  a 
sphere  of  radius  R  is  not  W,  but  slightly  less  than  V/, 


due  to  energy  losses 
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in  shock  wave  propagation. 

The  pressure  and  energy  at  a  point  (R,  0  )  has  been  evaluated  on  the 
basis  of  the  acoustic  approximation.  Now,  the  results  will  be  modified 
to  take  into  account  the  non-linear  nature  of  the  propagation.  This  pro¬ 
cedure,  while  not  being  rigorous  in  the  mathematical  sense,  does  give  a 
clear  picture  of  the  physical  phenomena,  as  well  as  numerical  results 
which  agree  quite  well  with  experimental  data. 


12.  Losses  Due  to  Shock  Wave  Propagation 


A  large  amount  of  work,  both  experimental  and  theoretical,  was  done 

prior  to  and  during  World  War  II  in  the  field  of  underwater  explosions 
20 

.  Several  theories  were  set  forth  for  the  propagation  of  the  shock  wave 

21  22 

formed  by  the  detonation  of  an  underwater  explosive  ’  "  . 


These  theories  may  be  applied  to  the  shock  wave  produced  by  an 

underwater  spark  discharge  because  of  the  "principle  of  similarity.  " 

The  principle  of  similarity  states  that  if  the  linear  size  of  a  charge  is 

changed  by  a  factor  k,  the  pressure  conditions  will  remain  unchanged  if 

new  distance  and  time  scales  k  times  as  great  as  the  Oi*iginal  ones  are 
23 

used  .  In  the  formulation  of  this  principle,  all  forces  not  scaling 
geometrically  have  been  neglected. 


For  spherical  symmetry,  the  general  characteristics  of  shock  wave 

propagation  are  such  that  if  the  shock  pressure  is  assumed  to  have  the 
f  I A  ’ 

form  pe  '  at  the  charge,  the  peak  pressure  will  decay  at  a  rate  more 

r  t 

rapidly  than  1/r  and  be  onlv  a.  function  of  R/a  ,  where  a  is  the  initial 

o  o 

charge  radius. 

I 

The  time  constant  ^  will  increase  with  R/a  ^  representing  a  broaden¬ 
ing  of  the  wave  as  it  propagates  outward. 

The  case  of  cylindrical  symmetry  is  considerably  more  complex  to 
analyze,  and  only  approximate  theories  have  been  developed.  If  we  re¬ 
strict  the  discussion  to  large  R/a^,  where  a^  is  the  initial  cylinder  radius, 
we  may  obtain  approximate  numerical  results. 

At  large  R/a^,  the  cylindrical  symmetry  should  become  unimportant 
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In  (64),  /Q  is  the  density  corresponding  to  the  pressure  p^  which  may 


be  determined  from  the  equation  of  state  for  water 


25 


(65) 


C  ‘ 
o  o 


It  is  interesting  to  note  that  the  asymptotic  time  constant  Q  does  not 
depend  directly  upon  the  time  constant  of  the  external  pressure  at  the 
plasma  water  boundary.  In  fact,  as  we  will  see,  ^  will  in  general  be 
many  times  greater  than  'T  . 


It  is  evident  that  the  peak  pressure  falloff  in  excess  of  acoustic  and 
the  time  spread  of  the  shock  wave  will  result  in  a  dependence  of  the  energy 
flux  upon  distance.  It  is  found  for  underwater  explosions  that  the  energy 
flux  is  not  as  dependent  upon  range  as  the  peak  pressure  since  it  re¬ 
presents  an  integration  of  the  entire  pressure  time  curve  and  the  effects 
of  large  peak  pressures  are  smoothed  out. 
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IV.  THE  ELECTRICAL  PROPERTIES  OF  THE  UNDERWATER 

SPARK  DISCHARGE 


1 .  An  Approximate  Equivalent  Circuit 

An  approximate  equivalent  electrical  circuit  for  the  underwater  spark 
discharge  may  be  obtained.  If  the  variations  of  plasma  conductivity  CT 
with  time  are  neglected,  the  electrical  resistance  of  the  plasma  is, 


(66)  R  R  /x^ 

g  o 


where. 


(67) 


L 


The  variation  of  x  with  t  is  given  by  (23).  If  R^  is  sufficiently  large 
compared  to  the  series  loss  resistance  of  the  discharge  circuit  and  the 
reactance  of  the  series  loss  inductance,  the  equivalent  circuit  of  figure  11 
may  be  drawn. 


t=0 


Figure  11.  Underwater  Spark  Discharge  Approximate 
Equivalent  Circuit 
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The  storage  capacitance  C  is  assumed  to  be  charged  to  a  potential  V, 
and  at  t  =  0,  breakdown  occurs.  The  following  differential  equation  may 
be  written; 


(68) 


Where  the  initial  condition,  q  (o)  =  CV  must  be  satisfied.  The  variables 
are  easily  separated  giving; 

(69)  dq  ,  ,  _  (  dt 


g 


Transforming  to  the  x  variable  from  (21)  and  using  (66); 


X 


(70) 


2  '2' 

logq=-^ 


IT 


dx 


After  performing  the  integration  ,  and  using  the  initial  condition,  we 
obtain  the  following  solutions: 


(71) 


V  =  Vexpf 


l^x^  (x^-1)^^^  -  3/4  X  (x^  -1)^^^  -3/6  X 


(x‘ 


,.1/2^3  , -1  ll 

- 1)  cosh  X  I  r 


and, 


(72)  i  =  V/R  x^  exp 


2  r  -2r  r  3  ,  2  ,^l/2  3  ,2  ,,3/2 

o  —  1  X  (x  -1)  x(x  -1) 


I 


u 


-3/8x(x  -I)  +3/8  cosh  ^ 
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Using  (24),  for  X  1.2,  t/T;^l.  3  these  expressions  become 


(73)  V  V  exp  (t/T:)2j 

and. 


(74) 


V/R  (t/T  )  exp 


[- 


'r/2R  c 

o 


(t/T) 


■] 


Equations  (73)  and  (74)  may  be  used  to  solve  for  the  instant  of  maxi¬ 
mum  current,  the  current  at  maximum,  the  voltage  at  the  instant  of  maxi¬ 
mum  current,  and  the  plasma  resistance.  The  following  relationships  are 
obtained: 


(75) 


/  t  Ir  cr 

m  g  o 


V  .  6  V 
in 


(a) 


(b) 


(c) 


C 


(d) 


By  referring  to  figure  3,  it  is  evident  that  the  form  of  (73)  and  (74) 
is  matched  to  the  experimental  current  and  voltage  waveform.s. 


2.  Pinch  Pressure 


The  magnetic  field  set  up  by  the  current  flowing  through  the  spark 
channel  produces  average  forces  on  the  charged  particles  tending  to  force 
them  together.  The  net  reduction  of  pressure  on  the  constraining  envelope 

2, 7 

produced  by  this  effect  is  known  as  the  pinch  pressure.  * 

If  we  assume  the  total  current  to  be  uniformly  distributed  over  the 
cross-sectional  area  of  the  plasma,  the  pinch  pressure  is, 
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where  (10)  henry  s/meter 


This  parabolic  pressure  distribution  when  integrated  over  the  channel 
cross-section,  gives  an  average  pinch  pressure  of; 


_  Mi;  ^  _  i^(lO)  ^  newtons 
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meter' 


If  t/'t  >  1.  3,  the  approximations  (25)  and  (74)  may  be  used  in  (77)  to 
obtain  the  pinch  pressure  as  a  function  of  time.  * 
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The  maximum,  pinch  pressure  occurs  at  a  time; 
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The  maximum  value  of  the  average  pinch  pressure  is; 


(80) 


-  _  .  55  (1U)~^  V‘ 

’’•’“"‘'max. 


newtons 
meter  ^ 


It  has  been  found  that  the  magnitudes  of  the  underwater  spark  discharge 
pinch  pressures  are  sufficiently  small  so  as  to  have  a  negligible  effect 
upon  the  observable  properties  of  the  underwater  spark  discharge. 


V,  COMPARISON  WITH  EXPERIMENT 


A  comparison  of  the  theoretical  results  with  experimental  data  will 
be  made  for  an  underwater  spark  discharge  having  the  parameters, 

V  =  30  kilovolts,  C  =  .  1  ^f,  L  =  4  cm,  E  =  45  joules. 
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These  pai  aiaetej.  b  were  selected  since  Frolov  and  Roi*'  provide 
measured  current,  voltage  and  pressure  waveforms  for  these  values. 

Current  and  Voltage  Waveforms 

When  the  parameter  R^^t^  in  equations  (75)  is  known,  the  theory  pre¬ 
dicts  the  expected  current  and  voltage  waveforms.  However,  Rq'^  must 
be  determined  from  one  of  the  experimental  values.  Since  the  time  be¬ 
tween  breakdown  and  current  maximum  appears  to  be  insensitive  to  minor 
variations,  the  experimental  value  of  t^  =  1.  1  yusec  will  be  used. 

Equation  (75a)  then  allows  us  to  calculate  =  12.  1  (10)  This  value 

has  been  used  in  (73)  and  (74)  to  calculate  the  theoretical  current  and 
voltage  waveforms  which  are  plotted  in  figure  12  along  with  the  experi¬ 
mental  curves.  The  differences  may  be  resolved  qualitatively. 
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Martin  determines  the  variation  of  plasma  conductivity  during  the 
discharge.  The  conductivity  is  found  to  increase  to  a  maximum  near  the 
instant  of  maximum  current  and  then  decrease.  If  the  plasma  conduct¬ 
ivity  behaved  in  the  same  manner  in  this  case,  it  would  explain  the  fact 
that  the  experimental  current  becomes  greater  than  the  theoretical  value 
near  the  current  maximum,  and  then  drops  more  rapidly  than  the  theoreti¬ 
cal  curve. 


In  order  to  evaluate  the  theoretical  results  for  the  acoustic  properties. 


a  value  of  X  is  necessary.  An  approximate  estimate  may  be  obtained  from 
the  experimental  current  and  voltage  values  after  breakdown,  which  are 
V 25  kilovolts,  I  ^370  amperes.  Dividing  these  values,  we  obtain  a 
value  of  =  68  ohms.  Using  the  previously  determined  value  of 

=  12.  1  (10)  a  value  of 't*  =  .  IS^sec  is  obtained. 

The  approximate  equations,  (73)  and  (74),  should  be  valid  for  x  ^  1.  2, 
t  ^  1.  3  ,  or  t  ^  .  235  ^sec,  which  covers  the  greater  part  of  the  current 

and  voltage  curves. 

Acoustic  Properties 

A  numerical  evaluation  of  the  theoretical  results  for  the  acoustic 
properties  of  the  underwater  spark  discharge  must  be  based  upon  esti¬ 
mated  values  of  ^  and  Y.  Since  the  electrode  voltage  drops  at  the  instant 
of  breakdown  from  30  kilovolts  to  about  20  kilovolts  at  current  maxim.um 
approximately  a  microsecond  later,  it  appears  reasonable  to  pick  the 
energy  represented,  25  joules,  as  a  value  of  Y.  This  value  of  Y,  together 
with  the  estimated  value  of  T  permit  (24)  and  (26)  to  be  solved  simultan¬ 
eously  for  p  and  a  .  The  calculated  values  are. 


p^  =  12,  200  atmospheres 

a  =  . 4  mm 
o 


These  values  may  be  used  to  check  the  assumption  of  incompressibil¬ 
ity  which  is  valid  provided  the  wall  velocity  is  much  less  than  the  speed 
of  sound.  The  calculated  niaximum  wall  velocity  is  from  (34); 
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555  meters / sec 


The  velocity  of  sound  is  1500  meters/sec,  so  that  the  maximum  wall 
velocity  is  about  1/3  that  of  sound.  If  it  is  remembered  that  this  maximum 
velocity  exists  over  a  short  part  of  the  expansion  of  the  plasma,  the  as¬ 
sumption  of  incompressibility  appears  reasonable. 


The  external  pressure  at  the  plasma  wall  calculated  from  (39)  is 
shown  in  figure  13.  The  fact  that  the  time  scale  is  much  less  than  that  of 
the  current  and  voltage  waveforms  is  evident. 


The  peak  pressure  which  should  be  measured  at  a  distance  R  from 
the  plasma,  in  a  direction  normal  to  the  channel  may  be  calculated  by 
means  of  the  shock  wave  relationships.  Using  (63); 
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1.1  (10)  (12,500)  =  13.  8  atmospheres 


Experimentally,  a  peak  pressure  of  9  atmospheres  is  measured  at 
that  distance  by  a  piezoelectric  probe  with  a  frequency  response  specified 
as  uniform  up  to  800  kc.  The  difference  between  the  calculated  and  ex¬ 
perimental  values  may  be  due  in  part  to  the  frequency  response  of  the 
probe. 


The  approximate  value  of  the  time  constant  at  this  distance  may  b( 
calculated  from  (64); 
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0  Cii  4.  45  ^sec 

The  measured  value  is  approximately  6^sec. 

Since  the  time  constant  at  the  plasma-water  boundary  was  calculated 
as  .  36  psec,  the  time  scale  appears  to  have  increased  by  about  a  factor 
of  12. 

Using  the  theory  for  the  directional  properties  of  the  pressure  wave¬ 
form,  the  amplitude  and  duration  of  the  waveform  in  a  direction  parallel 
to  the  spark  channel  may  be  calculated.  In  this  direction, 

_2 

=  L/c  sin  3  =  — ' — 5 —  =  27.  6  usee 
1.  5  (10)^ 

The  calcxilated  time  constant  at  0“  may  be  used  to  obtain  oC  "Lie, 


oc  L/c  = 


4  (10) 


-2 


=  5.9 


4.45  (10)  ^  1.  5  (10)' 


We  may  obtain  the  peak  pressure  at  90°  from  (49): 
P„  r  -  e-  ‘^'<=3 
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=  — g  ^  =  2.  34  atmospheres. 


A  value  of  2  atmospheres  is  measured  experimentally. 


The  experimental  and  theoretical  waveforms  for  directions  normal 
ai->d  parallel  to  the  underwater  spark  discharge  plasma  are  shown  in 

figure 


14. 
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EXTERNAL  PRESSURE 
AT  ONE  METER 


Fit?  u  re  14  Theor 


rheoretical  and  Experimental  Pressure -Time  Curves 


The  agreement  between  the  calculated  and  measured  values  is  better 
in  the  parallel  direction  since  the  probe  is  not  required  to  pass  a  steep- 
fronted  pressure  waveform,  as  it  is  in  the  normal  direction. 


Pinch  Pressure 


The  maximum  average  pinch  pressure  may  be  calculated  from  (80); 
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The  calculated  pinch  pressure  of  4.  1  atmospheres  is  sufficiently  small 
so  that  the  error  introduced  by  neglecting  the  effects  of  the  pinch  pressure 
upon  the  equation  of  motion  are  negligible. 
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Martin  calculated  the  pinch  pressure  for  an  underwater  spark  dis¬ 
charge  having  the  parameters,  V  =  25  kilovolts,  C  =  5.  8  ^f,  L  =  1.  5  cm. 

He  found  the  maximum  average  pinch  pressure  to  be  735  atmospheres  for  a 
peak  current  of  85,  000  amperes.  It  may  be  concluded  that  the  contribution 
of  the  pinch  pressure  to  the  kinetic  pressure  of  an  underwater  spark  dis¬ 
charge  is  not  significant. 

Electro-Acoustic  Efficiency 

The  experimental  energy  flux  directivity  factor  may  be  obtained  by 
calcuiating  the  energy  flux  radiated  in  a  direction  normal  to  the  discharge, 
and  relating  it  to  the  total  energy  radiated. 


K  the  waveform  measured  normal  to  the  discharge  plasma  were  the 
same  over  a  sphere  of  radius  R; 


4'nR^p'‘ 
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Frolov  and  Roi^  report  an  electro-acoustic  efficiency  of  32%  which 
they  calculate  by  integrating  pressure-time  curves  measured  over  a 
spherical  segment  of  1  meter  radius.  On  the  basis  of  the  definition  of  the 
energy  flux  directivity  factor  ^  from  (58); 


Therefore,  the  total  electro-acoustic  efficiency  should  be; 


The  measured  efficiency  should  be  lower  since  in  the  derivation  of 
(60)  shock  wave  energy  losses  were  neglected. 


The  good  agreement  between  the  calculated  and  measured  acoustic 


parameters  appears  to  justify  the  value  of  Y  -  25  joules,  which  was 
selected  as  the  energy  available  for  the  motion  of  the  spark  discharge 
plasma. 


VI.  SUMMARY 


A  theoretical  investigation  of  the  electro-acoustic  properties  of  the 
underwater  spark  discharge  has  been  presented.  Equations  for  the  con¬ 
servation  of  the  mass  and  momentum  of  an  incompressible  fluid  have  been 
used  to  derive  approximate  expressions  for  the  kinetic  and  potential 
energies  of  flow  outside  the  plasma- water  boundary.  Setting  the  sum  of 
these  energies  equal  to  a  constant  resulted  in  an  equation  of  motion 
involving  the  plasma  radius,  wall  velocity,  and  external  pressure. 

Assuming  adiabatic  expansion  of  the  plasma,  a  relationship  between 
the  external  pressure  and  plasma  volume  was  obtained.  This  relation¬ 
ship  permitted  an  explicit  solution  of  the  equation  of  motion  for  the  plasma 
radius  as  a  function  of  time.  The  solution  obtaijiied  v/as  expressed  in  terms 
of  the  ratio  of  the  plasma  radius,  a,  to  its  initial  radius,  a^,  and  the 
ratio  of  the  time,  t,  to  a  characteristic  time,'^,  which  was  found  to  be 
a  function  of  the  initial  values  of  plasma  radius,  and  external  pressure. 

For  a/a^  1  the  calculated  radius  becomes  almost  directly  proportional 

to  the  square  root  of  the  time,  a  result  which  is  in  agreement  with  ex¬ 
perimental  observations. 

The  calculated  plasma  radius  as  a  function  of  time  and  the  adiabatic 
pressure- volume  relationship  enabled  the  external  pressure  at  the  plasma 
wall  to  be  determined.  The  external  pressure  was  found  to  be  approxi¬ 
mately  an  exponential  function  of  time  with  a  time  constant  of  2 

An  exponential  approximation  for  the  external  wail  pressure  was  used 
in  the  derivation  of  the  directional  distribution  of  the  radiated  pressure 
pniH'-,  i''rom  the  agr eement  between  the  experimental  and  calculated 


pressure-time  curves  in  directions  normal  to  and  parallel  to  the  spark 
channel  axis,  one  may  conclude  that  to  a  good  approximation  the  under¬ 
water  spark  discharge  behaves  as  a  uniform  linear  acoustic  source 
radiating  an  exponential  waveform. 

Finite  amplitude  effects  have  been  discussed  through  application  of 
the  shock  wave  propagation  theory  developed  for  an  underwater  explosion. 
Major  effects  such  as  peak  pressure  falloff  in  excess  of  acoustic  and 
broadening  of  the  waveform  as  it  propagates  outward  have  been  discussed. 

By  neglecting  the  variation  of  the  plasma  conductivity  with  time,  and 
utilizing  the  equation  of  motion  for  variation  of  the  spark  channel  area, 
an  expression  for  the  plasma  resistance  has  been  derived.  This  resist¬ 
ance,  together  with  a  simple  electrical  equivalent  circuit  has  been  used 
to  obtain  theoretical  discharge  current  and  voltage  waveforms. 

These  results  show  that  measurable  quantities  associated  with  the 
experimental  current  and  voltage  waveforms  may  be  used  to  estimate 
parameters,  such  as  the  initial  energy  available  for  plasma  motion,  Y,  and 
the  characteristic  time'll  .  Once  they  are  known,  these  quantities  allow 
the  prediction  of  other  electro-acoustic  properties  of  an  underwater  spark 
discharge.  Quantities  such  as  initial  peak  wall  pressure  and  channel 
radius,  which  are  not  easily  obtainable  by  experimental  methods,  may  be 
calculated  by  utilizing  a  set  of  simultaneous  equations. 

Through  the  above  techniques,  the  initial  wall  pressure  and  plasma 
radius  for  a  30  kv,  .  1  ^fd,  45  joule  underwater  spark  discharge  have  been 
calculated.  The  values  obtained  for  p^  and  a^  were  respectively,  12,  200 
atmospheres,  and  .4  mm.  Using  these  values,  the  properties  of  the 
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radiated  pressure  pulse  have  been  determined  and  found  to  be  in  good 
agreement  with  the  experimental  results. 


VII.  PRACTICAL  APPLICATIONS 


A  number  of  areas  exist  in  which  the  underwater  spark  discharge  may 
be  of  great  practical  utility. 

At  the  present  time,  several  large  organizations  are  investigating  the 

application  of  the  underwater  spark  discharge  to  metal  forming.  A  n\imber 

of  difficult  forming  jobs  have  been  performed  quite  rapidly  by  placing  the 
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piece  to  be  formed  between  a  die  and  an  underwater  spark  discharge. 

Another  application  for  the  underwater  spark  discharge  is  a  sonic 
source  for  underwater  echo  sounding.  Specific  applications  are  classi¬ 
fied  information  and  cannot  be  mentioned;  however,  the  underwater  spark 

3 1 

discharge  shows  great  promise  in  this  area. 

The  underwater  spark  discharge  might  be  used  for  an  investigation  of 

the  elastic  constants  of  materials  under  high  transient  stress.  Through 

measurements  of  refracted  and  unrefracted  shock  fronts  at  a  water-metal 

boundary,  the  velocity  of  the  shock  wave  in  the  material  could  be  calculated, 

32  33 

and  the  elasticity  of  the  material  obtained.  ’ 


Vm.  SUGGESTIONS  FOR  FUTURE  STUDY 


There  is  a  great  amount  of  work  which  remains  to  be  done  before  the 
underwater  spark  discharge  will  be  as  well  understood  as  some  of  the 
other  types  of  gaseous  discharges. 

A  detailed  analysis  of  the  breakdown  of  a  water  gap  and  the  initial 
stages  of  plasma  formation,  from  both  a  theoretical  and  experimental 
point  of  view,  would  be  valuable.  Such  a  study  might  be  similar  to  the 
one  performed  by  Martin\  except  that  a  number  of  different  parameter 
combinations  would  be  employed.  This  investigation  should  include  a  de¬ 
tailed  determination  of  the  electrical  characteristics  of  discharges  oc- 
curing  in  water  of  controlled  temperature  and  salinity.  In  particular, 
the  variation  of  plasma  conductivity  with  time  should  be  measured.  This 
study  would  be  similar  to  the  one  performed  by  Martin,  except  that  a 
number  of  different  parameter  combinations  might  be  employed. 

In  addition,  an  experimental  investigation  of  the  pressures  which 
exist  in  the  vicinity  of  the  underwater  spark  discharge  channel  would  be 
of  interest.  However,  the  problems  connected  with  such  measurements 
are  formidable.  If  a  pressure  sensitive  device  were  used,  it  would  have 
to  be  small  enough  not  to  disturb  the  field;  have  a  bandw'idth  of  the  order 
of  several  megacycles;  be  insensitive  to  sti'ay  pickup  from  extremely 
strong  electromagnetic  fields;  and,  finally,  have  sufficient  strength  to 
withstand  the  tremendous  pressures  involved.  Needless  to  say,  such  a 
device  may  be  hard  to  come  by.  Other  techniques  might  be  useful.  For 
example,  a  measurement  of  the  curvature  of  a  light  ray  passing  through 
a  high  pressure  region  would  be  an  indication  of  the  pressures  involved. 


The  accuracy  of  such  a  measvirement  is  limited  by  a  knowledge  of  the  vari¬ 
ation  of  the  refractive  index  of  water  with  pressure  and  suitably  accurate 

34 

measurements  may  not  be  available. 

A  related  investigation  might  consist  of  a  determination  of  electro¬ 
acoustic  efficiency  as  a  function  of  voltage,  energy,  and  gap  separation. 

This  study  would  be  particularly  valuable  for  voltages  greater  than  30 

2 

kilovolts  and  gaps  longer  than  5  cm,  since  Frolov  and  Roi  stopped  their 
investigation  at  these  values. 
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